Mon. Not. R. Astron. Soc. 000,[T]{8](2009) Printed 7 May 2009 (MN MbX style file v2.2) 



Nonresonant Grain Acceleration in MHD Turbulence 
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ABSTRACT 

We discuss a new type of dust acceleration mechanism that acts in a turbulent magnetized 
medium. The magnetohydrodynamic (MHD) turbulence can accelerate grains through res- 
onant as well as nonresonant interactions. We show that the magnetic compression provides 
higher velocities for super- Alfvenic turbulence and can accelerate an extended range of grains 
in warm media compared to gyroresonance. While fast modes dominate the acceleration for 
the large grains, slow modes can be important for sub-micron grains. We provide comprehen- 
sive discussion of all the possible grain acceleration mechanisms in interstellar medium. We 
show that supersonic velocities are attainable for Galactic dust grains. We discuss the conse- 
quence of the acceleration. The implications for extinction curve, grain alignment, chemical 
abundance, etc, are provided. 

Key words: dust, extinction-ISM: particle acceleration-kinematics and dynamics-magnetic 
fields. 



1 INTRODUCTION 

Dust is an important constituent that is essential for heating and 
cooling of the interstellar medium (ISM). It interferes with observa- 
tions in the optical range, but provides an insight to star-formation 
activity through far-infrared radiation. It also enables molecular hy- 
drogen formation and trac es the magnetic field via emission an d 
extinction polarimetry (see lHildebrand etal, 2000; Laz arianl2007t) . 
The basic properties of dust (optical, alignment etc.) strongly de- 
pend on its size distribution. The latter evolves as the result of grain 
collisions, whose frequency and consequences (coagulation, cra- 
tering, shattering, and vap orization) depend on grain relative ve- 
lociti es (see discussions in iDrainell 19851; lYan 1 Lazarian. & Drainel 
|2004 hereafter YLD04; lfflrashita & Yanll2009l) . 

All these problems require our understanding of grain mo- 
tions in the interstellar mediu m. The interstellar medium is 
magnetized and turbulent (see lArons & Max! Il975t) . Although 
turbulence has be e n invoked by a number of authors (see 
Kusaka et al.l Il970l: IVoelk et al.l ll98Ct iDraind 1 19851 ; lOssenkod 
1993; Weid enschilling & Ruzmai kina 1994) to provide substantial 



grain relative motions, the turbulence they discussed was not mag- 
netized. Dust grains are charged, and their interactions with MHD 
turbule nce is different from t he hydrodynamic case. 

In lLazarian & Yanl j2002l) and Yan & Lazarian (2003, hence- 
forth YL03), we considered a cceleration of charged dust by the 
MHD turbulence. In particular, Lazarian & Yan (2002) applied the 
theory of Alfven ic turbulence (jG oldreich & Srid harlll995L hence- 
forth GS95, see ICho. Lazarian. & Vishniad 120031 for a review ) 
to grain acceleration through gaseous drag; YL03 identified a 
new mechanism for grain acceleration through gyroresonance. In 
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YLD04, we showed that the acceleration by turbulence domi- 
nate over other processes and can provi de supersonic speed to the 
grains. On the other hand, it is know n (Spitzer [19761 : IShull|[l97^: 
iMcKee et aill 19871 lEflison. Drurv. & Meverlll997^ that grains can 
also be efficiently accelerated through betatron process by the com- 
pressions at shocks. In this paper, we shall consider the betatron 
acceleration of dust by both incompressible and compressible tur- 
bulence. The intent is to provide a comprehensive view of current 
understanding of grain dynamics in MHD turbulence. 

To describe the turbulence statistics we use the ana- 
lytical fits to the statist i cs of Alfvenic modes obtained in 
ICho. Lazarian, & Vishniacl d2002h and compressible modes ob- 
tained in ICho & LazariarTfeOOl hereafter CL02), which is con- 
sistent with the spectra of turbulence velocity from observations 
(Lazarian 2008). In §2, we provide a brief description of the rel- 
evant properties of MHD turbulence. In §3, we study acceleration 
by large scale turbulence, where acceleration in both collisionless 
and collisional regimes are discussed. In §4, a review of alternative 
acceleration mechanisms is presented. In §5, we apply the results to 
interstellar grains and give a comparison of different mechanisms. 
Implications of our results are included in the discussion of in §6, 
while the summary is provided in §7. 



2 MHD TURBULENCE AND ITS STATISTICAL 
PROPERTIES 

Turbulent acceleration may be viewed as the acceleration by a spec- 
trum of MHD waves that can be decomposed into incompressible 
Alfvenic, and compressible fast and slow modes (see CL02). Sepa- 
ration of MHD perturbations in compressibl e media into fast, slow 
and A lfven modes is discussed in GS95, Lithwick & Goldreich 
< l200ll) . and CL02. The actual decomposition of MHD turbulence 
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into Alfven, slow and fast modes was performed in CL02, who also 
quantified the intensity of the interaction between different modes 
(see below). Unlike hydrodynamic turbulence, Alfvenic turbulence 
is anisotropic, with eddies elongated along the magnetic field. This 
happens because it is easier to mix the magnetic field lines per- 
pendicular to the direction of the magnetic field rather than to bend 
them. As eddies mix the magnetic field lines at the rate kj_v k > where 
kj_ is the perpendicular component of wavenumber k, v# is the mix- 
ing velocity at this scale, the magnetic perturbations (waves) prop- 
agate along the magnetic field lines at the rate IcuVa, where hi is the 
parallel component of the wave number k, va is the Alfven speed. 
The Alfvenic turbulence is described by GS95 model. The corner 
stone of the GS95 model is a critical balance between these rates, 
i.e., fcx v A- ~ ku Kli which may be also viewed as coupling of eddies 
perpendicular to the magnetic field and wave-like motions parallel 

to the magnetic field. From these arguments, the scale dependent 

2/3 

anisotropy k\\ oc £y and a Kolmogorov-like spectrum for the per- 
pendicular motions can be obtained. It was conjectured 
in Lithwick & Goldreich (2001) that the GS95 scaling should be 
approximately true for Alfven and slow modes in moderately com- 
pressible plasma. For magnetically dominated, the so-called low (3 
plasma, CL02 showed that the coupling of Alfven and compress- 
ible modes is weak and that the Alfven and slow modes follow 
the GS95 spe ctrum. This is consistent with the analysis of obser- 
vational data (Lazarian 2008) as well as with the electron density 
statistics (see [Armstrong, R ickett, & Sp angler 1995). According to 
CL02, fast modes evolves separately and are isotropic. 

Depending on whether the turbulence velocity on the energy 
injection scale L,„y, hV is larger or smaller than the Alfven speed, 
turbulence can be either super- Alfvenic (Ma = (SV/va) > 1) or 
sub- Alfvenic (Ma < 1). Turbulence is considered super- Alfvenic 
in intracluster medium. Although less certain, the interstellar turbu - 
lence may also be in this regime (Beck 2001; Padoan et al. 2004). 
SuperAlfvenic turbulence evolves with hydrodynamic Kolmogorov 
cascade till the scale L = L;„y(8V /v^) 3 = Lj n jM^, where the turbu- 
lence velocity becomes equal to the Alfven speed vj^\. Below this 
scale, magnetic field becomes dynamically important and the MHD 
cascade described above sets in. 

We do not deal with the unbalanced turbulence (see 
iLithwick et alj200llBeresnvak & Lazariarj2008l ; lchandraij2008h 
in this paper, as the degree of imbalance in compressible in- 
terstellar is unclear. We consider all the turbulent energy is in- 
jected from large scales and do not include additional energy in- 
jection on small scales arisi ng from compressions of cosmic rays 
jLazarian & Bere snvak 200^). Note that the Alfvenic turbulence 
and slow modes are different from waves as they are short lived, 
only fast modes may resemble wave turbulence. 



3 ACCELERATION OF GRAINS BY LARGE SCALE 
TURBULENCE 

There exists an important analog between dynamics of charged 
grains and dynamics of cosmic rays (CRs) and we shall modify 



the existing machinery used for cosmic rays to describe charged 
grain dynamic^. 

The energy exchange between the particles and the waves in- 
volves both resonant and nonresonant interactions. Resonant inter- 
action was considered in YL03 and YLD04. Here we shall be inves- 
tigating the role of nonresonant interactions on grain dynamics. The 
role of the non-resonant acce leration of CR by large scale compres- 
sions were discussed earlier dBerger et al.l ll958 ; Kulsr ud & Ferraril 



19711; lPtuskmlll98 8). and further clarified for MHD turbulence in 
Cho & Lazarianl J2006T) . but there are substantial differences be- 



tween cosmic rays and grains. In particular, the gyroradii of grain 
are much larger and can get comparable to the correlation length of 
turbulence. Moreover, grains are moving much more slowly than 
the CRs and therefore scattering is not important. Mean free path 
l m f p of grains are thus much longer and there exists a more impor- 
tant collisionless interaction regime with / < l m f„, which we shall 
demonstrate is the dominant nonresonant interaction in the paper. 

Apart from acceleration, a grain is subjected to gaseous fric- 
tion. For the sake of simplicity we assume here that grains are mov- 
ing in respect to stationary gas and turbulence provides nothing but 
the electromagnetic fluctuations. The acceleration of grains aris- 
ing from the gas motion is considered separately (see LY02). We 
describe the stoc hastic acceleration of dust by the Br ownian mo- 
tion equation (see Kal mvkov. Coffey. & Waldronll 996): mdu/dt = 
— u/S + Y, where m is the grain mass, u is the grain speed, Y 
is the stoc hastic accelera t ion force, 5 = td rn „lm is the mobility 
coefficient (Dra jneH"l985l ; lEazarian & Yanl l2002h . The drag time 
due to collisions with atoms is essentially the time for collisions 
with the amount of gas with the mass equal to that of the grain, 
''drag = ( a P^/-/' J «)\/ It / (8m n fc#r), where a is the grain size, T is 
the temperature, n„ is the neutral density, p gJ - is the mass density 
of grain. We adopt p gr = 3.8gcm~ 3 for silicate grains. The ion- 
grain cross-section due to long-range Coulomb force is larger than 
the atom-grain cross-section. Therefore, in the presence of colli- 
sions with ions, the effe ctive drag time decreases by the factor 
llDraine & Salpeteilll979h . 
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where rin is the hydrogen density, x; is the abundance, relative to 
hydrogen, of ion with mass m,-, x = Ya x U and fz is the probability 
of the grai n being in cha rge state Z. When the grain velocity gets 
supersonic ( |Purcellll969t) . the gaseous drag time is given by t s d = 

t drag/(^-T^ u / c s + 0.75cs/u — c^/2u 2 + Cr/5a 3 ), where c$ is the 
sound speed. 

Multiply the above Brownian equation above by u and take 
the ensemble average, we obtain 



d <u z > 

i 

dt 



< £ >, 



(2) 



Following the approach similar to that in Metros^ J 198Ct) . we can 
get from Eq.Q the energy gain rate < £ > for the grain 



Although in general turbulent generation of magnetic field which will 
bring kinetic and magnetic energy to equipartition, the dynamo process 
takes time and super-Alfvenic turbulence can appear as a transient phe- 
nomenon. 



2 In what follows we assume that the time scale for grain charging is much 
shorter than the grain Larmor period. 
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<e> = ——^ — (up Dp} ~ — , 
4p z op m 



(3) 



where D p is the momentum diffusion coefficient, and shall be de- 
rived below for different circumstances. 



3.1 Acceleration in incompressible turbulence 

One important nonresonant process is betatron acceleration. It 
arises from the compression of the magnetized fluid, which induces 
nonzero curl of electric field V x E in the plane perpendicular to the 
local magnetic field. The same betatron process was studied earlier 
in the context of shock acceleration, where grains are accelerated 
at postshock due to the compression. Here we shall investigate the 
betatron acceleration in turbulence. We stress that the betatron pro- 
cess operates even in incompressible turbulence, where magnetic 
field can be compressed by pseudo-Alfven modes, which are the 
incompressible limit of slow modes. 

Let us start with incompressible turbulence. The aver- 
age electric field dur ing one gyro-orbit can be estimated as 
dCho & Lazarianll2006l) . 



eddy turnover time and the particle streaming time. This is because 
the scattering of grains are infrequent because of the low speed as 
pointed out in YL03. 



3.2 Acceleration in compressible turbulence 

In compressible medium, there are both slow modes and fast 
modes. For compressible modes, the above mechanism also oper- 
ates on the scale / < l m fp- The compression on large scale / > l m fp 
is much less efficient in accelerating grains for the following rea- 
son. Since usually l m fp 2> r g , I > l m f p indicates / ~S> r g . Lack of 
sufficient scattering the adiabatic invariant of grain is con- 

served in this case and the ch ange of grain velocity is only a factor 
of few JCho & Lazarianll2006h . 

The betatron acceleration by fluctuations on I < l m t p does 
not work for small grain, whose mean free path is less than the 
turbulence damping scale l c . Lack of scattering, the mean free 
path of dust is usually determined by the gaseous or plasma drag 
*mfp = u ^drag- 

Let us consider the betatron acceleration by the compressible 
modes on the small scales / < l m f p . The first term in Eq.10 is 
nonzero and needs to be accounted for. 
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— [ [Vx (vxB)l -ds. 
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— f (BV-v-Bo-Vv)-ds, 
lr e c J 



(4) 



where r g is the Larmor radius of grain, c is the light speed. The first 
term arises from the compression of the fluid, which is equal to zero 
in incompressible fluid. We approximate the total magnetic field B 
by the mean field Bo in the second term. We also used the fact that 
the MHD wave in incompressible medium is transverse, i.e., v _L k. 
Therefore we focus here on the second term. In this case, 



E = 



2%r s I ^ 



VxE)-ds, 



2%r g c 



[Vx (vxB)] -ds, 



k±Vk,±Bpr g 
2c 



(8) 



where Vy- j_ is the fluid velocity perpendicular to the magnetic field. 
Accordingly, the electrical force and the momentum diffusion co- 
efficient are 



dp e 
~dt ~ 



■ k±v k \pj_/2, 



(9) 



E ~ — BohtVfc \\r g /2c. 



(5) 



, where v k || is the fluid velocity parallel to the magnetic field. The 
rate of momentum change is then 



dpi 
dt 



■- qE ~ qBok\\v kj \\r g /2c = p±k\\v k> i\/2. 



(6) 



where q = Ze is the charge of grain and pj_ is the component of 
grain momentum perpendicular to the magnetic field. The accel- 
eration is canceled out by the adiabatic loss in the parallel direc- 
tion if the particle velocities are isotropic. Thus it only operates o n 
the scale less than the mean free path l m fp 

Lazarian 2006), 

from which we get 



At- 



'Pi(V*,||) 2 /4xAr ' 



6Lmax(v / t,i<) 



(7) 



where we adopted the GS95 scaling for the pseudo-Alfven modes, 
v k ~ V(k^L)~ 1 /^ and ku ~ £-^ 3 L~'/ 3 , V is the turbulence veloc- 
ity at the injection scale L of MHD turbulence. For pseudo Alfven 
modes, _L k and resides in the k B plane. Morever most en- 
ergy is concentrated in the modes with k _L B in the GS95 model. 
Therefore we used v k n ~ v k in the above equation. Unlike cosmic 
rays, At = l/[fei|max(vA,«)] is determined by the minimum of the 



A? ~ k 2 ± p 2 ± v 2 k x nmx(t wav ,t gc )/4, 



(10) 



where t waY is the wave period and t gc is the confinement time of the 
grains. 

For slow modes, 



Dp± 



kTL-^r?. V 2 



4max( 



'ph>u) v 2 k 



(11) 



where v p h = min(vyi,cs) is the phase speed of the slow modes. The 
velocity displacement of slow modes is in the direction of kii — 

6 sin 8 cos 8 k^ and _L k in the high (3 medium dCho & Lazarianl 
2006), where 6 is the wave pitch angle, i.e., the angle between 
the wave vector k and B, p = 2c|/v^ is the ratio of thermal pres- 
sure to magnetic pressure. Therefore, < v 2 . ± /v 2 >~< cos 2 > in 
high p medium and < v\ L /v\ p 2 < sin 2 6 cos 2 6 > /4 in low P 
medium, Because of the strong anisotropy (k± 3> fell), < sin 2 6 >~ 
1, and < cos 2 G >^ (k\\/k L ) 2 = (/t ± L) -2 / 3 . All the scales from 
lmfp to the truncation scale l c contribute to the acceleration. Since 
the contribution from the scales smaller than the gyroradius has 
been included in the treatment of gyroresonance in YL03, we only 
account for the eddies larger than gyroradii. Integrating Eg. ill It. we 
get 
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Figure 1. (a): the compression of magnetic field induce a nonzero curl of electric field, which accelerate grains as grains are orbiting around the magnetic 
field; (b): the direction of velocity perturbation and electric field induced by compressible modes. 
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(12) 



In low P medium, the slow modes cause marginal compres- 
sion of magnetic field with the velocity perturbations being nearly 
parallel to the magnetic field. That is why the acceleration by the 
slow modes in the low p regime is substantially reduced. 

For fast modes, the velocity perturbation is ± B in the low P 
case and || k in the high p case and therefore v? , /v? ~ 1 for p < 1 
and 2/3 for P > 1 since fast modes are isotropic. Unlike slow modes, 
the cascade of fast modes is slower, namely, the cascading time is 
longer than the wave period. Therefore, one needs to take into ac- 
count periodic motion of the fast modes. In case when particles 
are confined in the eddy longer than the wave period, the expan- 
sion ensuing the compression will take the energy back from the 
particles. Therefore, the acceleration only happens if the particles 
leave the eddies within one compressional period. Since grains ex- 
perience marginal scattering, it indicates that only super-Alfvenic 
grains (p < 1) or supersonic grains (P > 1) can be accelerated. 



4 ALTERNATIVE ACCELERATION MECHANISMS 
4.1 Gyroresonance 

Gyroresonance acceleration of charged grains by a spectrum of 
MHD waves decomposed into incompressible Alfvenic, and com- 
pressible fast and slow modes (see CL02) was first described in 
I Yan & Lazarianl d2003h . Gyroresonance occurs when the Doppler 
shifted frequency of the wave in the grain's guiding center rest 
nSl is a multiple of the grain gyrofre- 



frame CO. 



CO 



quency, un is the component of grain velocity parallel to the lo- 
cal magnetic field. For low speed grains, we only need to con- 
sider the resonance at n = 1. From the resonance condition, we see 
that gyroresonance requires perturbations on the Larmor scale and 
therefore only applicable to large grains whose gyroradii are larger 
than the dissipation scale of the turbulence l c . The gyroresonance 
changes both the direction and absolute value of the grain's mo- 
mentum (i.e., scatters and accelerates the grain). The ratio of the 
scattering and acceleration rates is approximately w 2 / v l' ^ or ' ow 
speed grains, therefore, acceleration is the main effect. 



dt 



At . 



tg C /4 ~ k 



- 2 ^L- 2 ^p 2 ± V 2 /{6u) 



(13) 

where we have used Vf, = (fcL)~^, i; = 1/3 for Kolmogorov turbu- 
lence and 1/4 for acoustic turbulence. Unless otherwise specified, 
we adopt i; = 1/4 to enable comparison with our earlier results 
(YL03, YLD04). Since it increases with the wave number k, the 
contribution on the smallest possible scale is dominant. Excluding 
the contributions from scales < r g , the corresponding scale is then 
max (r g ,l c ). Accordingly, 



2 1/2 



p z V 



9umax{r g J c y- 2 ^L 2 ^' 



(14) 



Comparing the above equations, we see that usually the accel- 
eration by fast modes is much more efficient than by slow modes. 
Low speed grains, however, can only be accelerated by the slow 
modes for the reason we discussed above. 



4.2 Transit time damping 

In the calculations in Sj3] we neglected the mirror force —p.VB = 
—mu 2 l S7B/(2B), which scatters the grains, but does not increase 
the energy. In the small amplitude limit SB <C B, particles should 
keep in phase with the wave /c||V| = CO in order to have effective 
interactions with the compressions of magnetic field. This pro- 
cess is often refe rred as transit time damping (TTD) in the lit- 
erature (see, e.g.. lAchterberd 1 198 lk ISchlickeiser & Milled 1 19981: 
lYan & Lazarianil2004l) . In this case, the mirror force redirect the 
energy gain in the perpendicular direction to the parallel direction 
without changing the total grain speed. In fact, our test calcula- 
tion of the TTD using quasilinear approach is consistent with the 
results in fj3] within order of unity. I therefore do not think that in 
terms of the acceleration TTD is different from the betatron process 
discussed in this paper. 
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4.3 Hydrodynamic drag 

In hydrodynamic turbulence, the grains motions are caused by the 
frictional interactions with the gas. On large scales grains are cou- 
pled with the ambient gas, and the slowing fluctuating gas mo- 
tions will only cause an overall advection of the grains with the 
gas (Draine 1985), which we are not interested. On small scales 
grains are decoupled. The largest velocity difference occurs on the 
largest scale where grains are still decoupled. Thus the character- 
istic velocity of a grain with respect to the gas corresponds to the 
velocity dispersion of the turbulence on the time scale frf rag . In the 
MHD case, the grain perpendicular motions are constraine d by the 
Larmor gyration unless t^^Q. < 1 dLazarian & Y an 2002). 

The grain motions get modified when the damping time scale 
of the turbulence t c is longer than either f^ rag or 1 /SI. In this case, 
a grain samples only a part of the eddy before gaining the velocity 
of the ambient gas. The largest shear appears on the damping scale, 
and thereby u ~ v c min(l/f2,^ rag )/x c , where v c is the velocity on 
damping time scale x c . 

4.4 Ponderomotive force 

A cross-fie ld acceleration by the p onderomotive force was sug- 
gested by IShukla & Stenflol §005) as an another independent 
mechanism in turbulence to provide grains with supersonic speeds. 
If this were true, it would indeed indicate a dramatic shift in the 
whole scenario, especially for small grains. We find however their 
formalism is erroneous. Let us start from the electric field that could 
arise in the varying magnetic field. Considering the momentum 
equation of the electron, where the electron inertia can be neglected 
in the low frequency MHD regime, 

n e eE = - — y e x B-y e T e Vn e , (15) 
c 

where n e , \ e , T e , y e are the density, velocity, temperature and adi- 
abatic index of electrons. Combining it with the quasi-neutrality 
condition n e = iij + Znj and the induction equation, 

4ne 

VxB= (niVi+ZndU — n e v e ), (16) 

c 

we obtain 



words, there is no such an addit i onal ponderomotive acceleration 
as discussed in Shukla & Stenflol J2005I) . 



4.5 Weak shocks 

Simulations indicate that the density structures in super- 
sonic MHD turbulence a re ass ociated with slow shocks 
jBeresnvak. Lazarian. & Choi 120051) . Therefore in supersonic 
turbulence, one also needs to take i nto account th e contribution 
of grain acceleration fr om shocks tepsteinl Il980l : iMcKee et al.l 
ll987l : lEUisonetal.lll997l) . The shocks in the supersonic turbulence 
require better quantitative description, nevertheless. 



4.6 Hi rocket and radiative pressure 

A different mechanism of driving grain motions is a residual im- 
balance in "rocket thrus t" between the opposite surfaces of a rotat- 
ing grain dPurcellll 1979b . This mechanism can provide grain rela- 
tive motions and preferentially move grains into molecular clouds . 
Three causes for the thrust were sue gested by IPurcelll dl979h : 
spatial variation of the accommodation coefficient for impinging 
atoms, photoelectric emission, and H2 formation. The latter was 
shown to be the strongest among the three. The uncompensated 
force in this case arises from the difference of the number of cat- 
alytic active sites for H2 formation on the opposite grain surfaces. 
The nascent H2 molecules leave the active sites with kinetic energy 
E and the grain experiences a push in the opposite directions. H2 
thrust potentially may be important for small grains. In the typical 
interstellar environment, small grains flip frequently according to 
recent study (Hoang & Lazarian 2009). As a result, the thrusts are 
smeared out. 

Dust grains exposed to anisotropic interst ellar radiation 
fields are subjected to forces as well as torques (PurcefH [l979l : 
Wein gartner & Draind l2001al) . These forces arise from photo- 
electric emission, photodesor ption as well as radiation pressure. 
Weing artner & Draind d2001al) considered the forces and calculated 
the drift velocity for grains of different sizes. The velocities resulted 
from these processes, however, are smaller than those arising from 
interactions with turbulence in the interstellar medium (YLD04). 



1 /VxB \ 

E = n,v; - Zn d u x B - y e T e Vn e (17) 

n e e \ 471 / 

where iij, v; are the density and velocity of ions. If we neglect the 
thermal term in the low P medium we consider, the electric field is 
perpendicular to the magnetic field and cannot accelerate the par- 
ticles longer than the Larmor period. The net effect of the electric 
field is to induce a drift in the direction E x B. The corresponding 
drift speed is 



ExB 



ck x SB 

4ltn e e 



(18) 



Since the acceleration is on the time scale of Larmor period, 



the wavenumber in the above equation is given by k : 



and ion 



velocity that should be adopted is the velocity of the ions (gas) on 
the Larmor scale. It turns out then the second term is much less 
than the first term in the above equation, and the drift is in fact 
the same proces s as th e gaseous drag we already accounted for in 
lLazarian & Yai] d2002l) and also the Fig.l in this paper. In other 



5 RESULTS 

To compare with earlier work (YLD04), we consider first the 
warm ionized medium (hereafter WIM) with the same parameters, 
T = 8000K, n e = 0.1cm" 3 , B = 3.35/jG (see tabld5}. Here, we 
consider grains in the range of 10" 7 cm< a < 10- cm, which carry 
most grains mass in ISM. We adopt the mean grain charge as given 
in YLD04, which was obtained by balancing the collisions with 
electr ons and collisions with ions as well as photoelectric emis- 
sion dWeingartner & D raine 2001b). We assume that for MHD tur- 
bulence the injection of energy happens at the scale L where the 
equipartition between magnetic and kinetic energies, i.e., V = V^, 
is reached. If we adopt a velocity dispersion 8V = 40km/s at the 
scale Z = 30pc, then turbulence in WIM is super-Alfvenic at this 
scale. If turbulence at this regime follows hydrodynamic cascade, 
then L is 3.75pc where the injection velocity is V = Va = 20km/s. 
The acceleration is dominated by the interactions with fast modes. 
Since fast modes are long lived and can only accelerate particles 
whose crossing time is shorter than the wave period, the betatron 
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acceleration only operates with fast-moving large grains. The re- 
sult for the WIM is demonstrated in Fig[2] where the outcome from 
different grain acceleration processes are compared. 

Next we consider a warm neutral medium (hereafter WNM), 
with T=6000K, n H = 0.3cm -3 , B = and the same tur- 

bulence injection velocity and scale as in WIM (see table 1 for 
other parameters). In partially ionized medium, a viscosity caused 
by neutrals results in decoupling on the characteristic time scale 
jLazarian & Yanll2002h tdamp ~ V^*AT 2 ~ (/«v„) _1 Ar 2 , where V„ 
is the kinetic viscosity, l„ is the neutral mean free path, v„ is the 
thermal velocity of neutrals. When its cascading rate T^T 1 equals 
to the damping rate t^ amp , turbulence is considered damped. This 
defines the truncation scale of the turbulence l c . = kV 2 /v p h for 

fast modes dCho et alj|2002h and Zt = fc j /3 L~'/ 3 V4 for Alfven and 
slow modes ( Gold reich & Sridhajf l995). Then for WNM, one gets 
l c ~ 8.3 x 10 15 cm for the fast modes and l c ~ 5.3 x 10 17 cm for 
Alfven and slow modes. This introduces a cutoff grain size for the 
acceleration by fast modes, for which the Larmor period is equal to 
the damping time scale. The result for grain acceleration is shown 
in Fig[2] As we see, the betatron acceleration operates for an ex- 
tended range of grains and can be more efficient compared to gy- 
roresonance. In particular, fast modes accelerate further the super- 
Alfvenic grains and slow modes dominant the acceleration for the 
smaller grains. 

In YL03, we proposed gyroresonance as one important accel- 
eration mechanism for grain dynamics. We showed that the mag- 
netostatic approximation for turbulence is invalid for grains, whose 
velocities are less or comparable to the velocity perturbations in 
turbulence. As a result, the resonance function in the quasi-linear 
theory is broadened from the resonance wavenumber kit res ~ 1/ r g 
to all k in the inertial range. To cut off the spurious resonance con- 
tribution from the large scale perturbations, we truncated the inte- 
gration for the momentum diffusion at the k« res and discarded all 
the interactions from the larger scales. A complementary study is 
conducted in this paper to include the acceleration from the large 
scales. Our results show that the nonresonant interaction is com- 
parable to the gyroresonance and can get even more efficient in 
super-Alfvenic turbulence and the strength of magnetic field. 

Our results show that only large grains can be effectively ac- 
celerated by the large scale compressions in turbulence. It is due to 
the fact that the scattering of grains is marginal. With the ballistic 
motions, the adiabatic invariant is conserved provided that the grain 
gyroradii are much smaller than the turbulence injection scale. The 
grain energy gain is thus limited. For large grains, their gyroradii 
are close to perturbation scale and there is no such a constraint. Our 
results are therefore self-consistent. 

How would the results vary as the parameters of the par- 
tially ionized medium vary? The relative importance of the gy- 
roresonance and betatron acceleration depends on the magnitude 
of the magnetic field as well as the intensity of turbulence. With 
the increase of magnetic field, the relative efficiency of betatron 
acceleration decreases compared to gyroresonance. On the other 
hand, the betatron acceleration increases faster with the intensity of 
turbulence than gyroresonance. More specifically, gyroresonance 
dominates in equipartition turbulence and sub- Alfvenic turbulence 
(Ma = SV /va < 1). In super- Alfvenic turbulence (M A = SV/va > 
1), betatron acceleration becomes important and works for ex- 
tended range of grains being less constrained by the turbulence 
damping. In Sj5] for instance, we adopt 8V = 40km/s at / = 40pc. 
If one takes 8V = 20km/s as in YLD04, the nonresonant and reso- 
nant acceleration rates are comparable. It is therefore sufficient to 



consider only the gyroresonance. Although turbulent generation of 
magnetic field tend to bring kinetic and magnetic energy to equipar- 
tition, the time re quired is long (a few crossing time at the injection 
scale according to lCho et alj2008T) . during which grains can be ac- 
celerated through the betatron process with the super- Alfvenic tur- 
bulence. Even if the velocity dispersions that we adopt in table. 1 
are excessive for the entire ISM, the local driving may bring them 
to such and even larger values. 



6 DISCUSSIONS 

This paper intends to encompass all possible situations where 
dust grains can be accelerated. As presented in $4] various 
processes have been proposed in the literature for the grain 
acceleration. In regard to the pla sma processes alone, there 
have been shock ac celeration (e.g., lEpsteirJ Il980l : iMcKee et al.l 
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1987), gyroresonancei Yan & Laz arian 2003), ponderomotive force 
( Shuk la & Stenfldl2005h . We intend to supplement the list by tak- 
ing into account nonresonant betatron acceleration in MHD turbu- 
lence. This proc ess was discu s sed in the context of co smic ray ac- 
celeration (e.g.. |Ptuski n 1988; C ho & Lazarianl l2006). but has not 
been considered for charged grains. The study here demonstrates 
that the betatron acceleration can be important in super-Alfvenic 
turbulence. 

The steady state for the strongly super-Alfvenic turbulence is 
not achievabl e due to dynamo generation. However, numerical sim- 
ulations (see ICho et alj|2008l) show that the steady state may not 
necessarily have the exact equipartition and kinetic energy may 
sometimes dominate. In fact, even in the situation when turbu- 
lence in terms of its mean velocity dispersion in the volume is 
sub- Alfvenic, density flu ctuations may induce super-Alfvenic tur- 
bulence. This is shown in lBurkhart et al. (2008). 

It has been shown that the composition of the galactic cos- 
mic ra y seems to be better correla ted with volatility of ele- 
ments I Ell ison. Drurv. & Mevedfl997l) . The more refractory ele- 
ments are systematically overabundant relative to the more volatile 
ones. This suggests that the material locked in grains must be 
accelerated more efficiently th an gas-phase ions (Epsteir] [T980l: 
Ellis on. Drurv, & MeveJ Il997l) . The stochastic acceleration of 
grains, in this case, can act as a preacceleration mechanism. The 
accelerated grains are more efficiently scattered by MHD perturba- 
tions, which is a key component for diffuse shock acceleration. 

Grains moving supersonically can efficie ntly vacuum-clean 
heavy elements as suggested by observations (Wakk er&Mathis! 
2000). Gr ains can also b e aligned if the grains get supersonic (see 
review byE azana nl2007l) . The supersonic grain motions will result 
in grain alignment with long axes perpendicular to the magnetic 
field. In fact, for irregular grains with nonzero degree of helicity, 
ev en subsonic motions can g et grains aligned as recently illustrated 
byE azarian & Hoand 120071) . The critical velocity is, however, still 
unclear. If taking into account mirror reflection, the energy gain of 
grains is not strictly in the perpendicular direction. Nevertheless, 
fast moving grains are aligned with the long axes perpendicular to 
the magnetic fiel d regardless of its directio n of motions according 
to recent studies jLazarian & Hoangll2007h. 

In our recent paper [tlirashita & Yan ( 2009]), we calculated 
the equilibrium grain size distribution based on our earlier result 
(YLD04) and obtained corresponding extinction curve in various 
phases of ISM. There we applied the gyroresonance, which is dom- 
inant in equipartition turbulence and subAlfvenic turbulence. We 
believe the bulk of interstellar medium has Ma < 1 because of the 
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T(K) 


rag (cm 3 ) 


n e (cm 3 ) 


Guv 


B(/jG) 


L(pc) 


V = va (km ■ ■ 




Z c (cm) 


WIM 8000 


0.1 


0.0991 


1 


3.35 


3.7 


20 






WNM 6000 


0.3 


0.03 


1 


5.8 


3.7 


20 




8.3 x 10 15 



Table 1. The parameters of idealized ISM phases and relevant damping. Among them, «// is the number density of H, n e is the number density of electron, 
Guv is the UV intensity scale factor, V is the injection velocity. L is the injection scale of fast modes, which is taken as the scale where the turbulence velocity 
is equal to the Alfven speed (see the text for details). The dominant damping mechanisms for fast modes are given with the corresponding damping timescale 
% c . WNM=warm neutral medium, WIM=warm ionized medium. 
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Figure 2. Grain velocity vs. size owing to different acceleration processes in WIM (left) and WNM (right). The result from betatron acceleration is represented 
by the dashed line. The '.' line represents the gyroresonance with fast modes. The dash-dot line is the result from hydro drag with fast modes. The dotted line 
refers to the threshold velocity for coagulation. Contributions from different processes are approximately additive as shown by the solid line. 



dynamo process, the basic conclusion there is valid therefore. In 
particular environments where turbulence gets super- Alfvenic, the 
process discussed in this paper should be accounted for. As I dis- 
cussed earlier, the betatron acceleration is most efficient for large 
grains (> a few tenths nm), so the effect is mainly on shattering. 
Depending on the intensity of turbulence and the strength of mag- 
netic field, the velocity can be increased by a factor of a few through 
the betatron acceleration, which could decrease the upper limit of 
the grain size in the warm media by a factor of few. In cold media, 
including cold neutral medium and the molecular clouds, the be- 
tatron acceleration is less efficient, for two reasons. First, the low 
levels of UV and low temperatures result in reduced grain charge. 
Secondly, because of the increased density, the frictional drag is 
increased. The betatron process is subdominant unless in highly 
super-Alfvenic turbulence, which could appear in a transient state. 



7 SUMMARY 

We calculate the acceleration of charged grains arising from be- 
tatron acceleration by MHD turbulence. We used the CL02's re- 
sults for the description of MHD turbulence, according to which 
Alfvenic turbulence and slow modes follows GS95 scaldings while 
fast modes are isotropic. We obtained estimates of grain velocities 
in interstellar medium. We showed that 

(i) MHD turbulence dominate grain acceleration in interstellar 
medium. 

(ii) The relative importance of betatron acceleration and gy- 
roresonance depends on the intensity of turbulence. While gyrores- 



onance dominates in sub- Alfvenic turbulence, betatron acceleration 
can be important in super-Alfvenic turbulence. 

(iii) Nonresonant acceleration with fast modes are important for 
large grains in warm media. Slow modes can accelerate sub-micron 
grains through the nonresonant interactions. 

(iv) Grain velocities depend on the local conditions of turbu- 
lence. 

(v) Dust gets supersonic via interactions with interstellar turbu- 
lence. The velocities obtained are sufficiently high to be important 
for shattering large grains and efficiently absorbing heavy elements 
from gas. The acceleration can also result in mechanical alignment 
of grains perpendicular to the magnetic field. 
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